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In presence of a strong external magnetic field the dynam- 
ics of solitons in a one-dimensional easy-plane Heisenberg an- 
tiferromagnet exhibits a number of peculiarities. Dynamics 
of internal soliton degrees of freedom is essentially quantum, 
and they are strongly coupled to the "translational" mode of 
soliton movement. These peculiarities lead to considerable 
changes in the response functions of the system which can be 
detected experimentally. 
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I. INTRODUCTION 

Study of soliton excitations in onc-dinicnsional (Id) 
systems have been attracting much attention, since the. 
pioneering work by Krumhansl and SchriefferJd Mikeskatl 
was the first one who calculated the contribution of soli- 
tons into observable quantities for the reaUstic system, 
viz. quasi-ld easy-plane ferromagnet CsNiFs in external 
magnetic field. His prediction of the soliton-induced cen- 
tral peak in thcpdynamical structure factor was confirmed 
experimentallyjj which gave rise to intensive theoretical 
and experimental studies. Since then, soliton signatures 
were reliably observed in static and dynamic character- 
istics of various Id magnets with both ferCPP- and antifer- 
romagnetic coupling, see recent reviews.clEl 

Easy-plane Id antiferromagnetic compound TMMC 
((CD3)4NMnCl3) in exterpal field is one of the most 
widely studied materialsElLfl (see also Ref. and refer- 
ences therein). Its magnetic ions have spin S — 5/2, 
and thus it is believed that the dynamics of TMMC 
can be well described on the basis of classical equa- 
tions of motion for the sublattice magnetization. It is 
well known that in presence of a strong magnetic field 
the equations of motion for antiferromagnet are rather 
complicated and cannot be reduced to the sine-Gordon 
(SG) model|-because of the importance of out-of-plane 
deviations .BEj Effective field-induced anisotropy com- 
petes with the easy-plane anisotropy, so that at certain 
critical field He — {HeHa)^^^ (where Hf, and Ha are 
the exchange and anisotropy fields, respectively) in-plane 
and out-of-plane magnon modes have equal excitation 
energies, and the SG approximation becomes obviously 
inadequate for H close to He- 

At H = He the "easy plane" effectively changes. There 
are, respectively, two types of static kink solutions,a 



one of which has lower energy at H < He and the 
other at if > He- Corresponding dynamic 
are much less understood. Different authorajll^ 
obtained contradictory results concerning their stabil- 
ity. It was shownHEjO that the effective plane of ro- 
tation of the sublattice magnetization in a kink depends 
on the kink velocity, and the two types of kinks trans- 
form into each other at some critical velocity Vc, Vc ^^ 
at H ^ He. 

The main feature of the problem, determining kink 
stability and dynamics in strong magnetic fields, is the 
preseacc-of out-of plane magnon mode localized at the 
kink.t£HlZrQ This internal mode is highly nonlinear at 
V -^ Vc and strongly coupled to the "translational" (zero) 
mode of the kink motionj3 which makes its analysis rather 
complicated. Moreover, it can be shown that because of 
the relatively small effective mass of solitons [the ratio 
of soliton and magnon masses is S for antiferromagnet, 
while for ferromagnets this ratio is proportional to the 
large parameter [He/HaY^^] the dynamics of internal 
soliton degrees of freedom in antiferromagnets is essen- 
tiallij-quantum even for "almost classical" values of S like 
5/20113 and exhibits subtle-eflJacts distinguishing between 
integer and half-integer S'.Ejcj Thus, to our opinion, the 
internal dynamics of solitons in TMMC in a strong-field 
regime can be highly nontrivial, and its consistent de- 
scription, as well as experimental observation, is a fasci- 
nating problem. 

In this paper we present a detailed theoretical study of 
the soliton dynamics in Id "almost-classical" spin-S* easy- 
plane antiferromagnet in strong external magnetic field 
close to the critical value He- We show that quantum 
effects (tunneling between two energetically equivalent 
kink states) lead to appearance of a new localized mode 
which is also coupled to the translational mode and thus 
depends on the soliton momentum P. At P = the tun- 
neling is suppressed for half-integer S due to the topolog- 
ical effects, but at nonzero P it becomes possible for any 
S. We calculate the contribution of the localized modes 
into the dynamic structure factor (DSF) and show that 
it leads to an additional peak concentrated at nonzero 
frequency, which can be detected by electron spin reso- 
nance (ESR) or inelastic neutron scattering (INS) tech- 
nique. Numerical estimations of the corresponding reso- 
nance frequencies and the peak properties for TMMC are 
given. To the knowledge of authors, at present the only 
known experiment probing internal soliton modes is the 
so-called "soliton magnetic resonance" Ej in the Ising-type 



quasi-ld antiferromagnet CsCoCla, and no such data ex- 
ist for the Heisenberg magnets. 

The paper is organized as follows: In Section O we 
introduce the model. In Sect. Ill we derive an effective 



model of internal kink dynamics, perform classical linear 
analysis of the spectrum of excitations against the back- 
ground of kinks and discuss the soliton stability problem, 
and in Sect. IV we perform quantization of the effective 
model and discuss the behavior of the energy levels of lo- 
calized modes. In Sect. ^ the soliton contribution to the 
response functions is considered, and numerical estima- 
tions of predicted effects for TMMC are given. Finally, 
Sec. VI contains concluding remarks. 



II. MODEL 

We consider a one-dimensional two-sublattice Heisen- 
berg antiferromagnet (AFM) of the TMMC type in ex- 
ternal magnetic field. This system is well described by 
the following Hamiltonian: 

'H = 2_^ \JSnSn+l 

n 

+ D,{S:,f - D^iSyf - g^iBHSn] , (1) 

where J > is the exchange constant, jjlb is the 
Bohr magneton, g is the Lande factor, Di and D2 are 
anisotropy constants, Di ^ D2 > 0, H is the magnetic 
field, and spins Sn are treated as classical vectors. In 
absence of the field {xy) is an easy plane, and y is the 
easiest axis in this plane (the case of "pure" easy-plane 
AFM corresponds to D2 = 0); the chain axis is along z. 
For TMMC the value of the spin is 5 = |, and the follow- 
ing values for the parameters of the classical Hamiltonian 
have been estabhshed£3 JS"^ = 85 K, DiS^ = 1.9 K, 
028"^ = 0.022 K, and g = 2.01. 

Long-wavelength dynamics of AFM can be described 
within the nonlinear sigma model approach. It is con- 
venient to introduce the antiferromagnetism vector (the 
sublattice magnetization) li — {821+1 — 82i)/'2S at each 
magnetic elementary cell. At low temperatures, when the 
magnetization rhi — {S2i+i + S2i)/'2S is small, I can be 
regarded as a unit vector, and after passing to the con- 
tinuum limit the following Lagrangian can be obtained: 

L = Ij^^ajdz |1(9JT - {djf ^Jl - ^ 



-(t-Hf+'^H-itxdtV) 



(2) 



— / dzl-iyixdtl) . 



Its derivation from different points of view the reader 
can find in Refs. p^-p5|; this approach is closely related 
to the commonly used Mikeska's formulationJj Here a is 



the lattice spacing, c — 2JSa/h is the limiting veloc- 
ity of spin waves, 7 = gfiB/h is the gyromagnetic ra- 
tio, xq = a[J/2{Di + 02)]^/'^ and Xb = a{J 1202^1'^ are 
the characteristic length scales. The last term in (||) is 
the so-called topological term which can be rewritten as 
271118 Q, where Q is the Pontryagin index (the winding 
number) of a given space-time configuration of the field 
I. This term is irrelevant in classical analysis, but is very 
important for the quantum theory since it distinguishes 
between integer and half-integer S. 

Within this approximation, the vector of magnetiza- 
tion TO is a "slave" variable and can be expressed through 
the antiferromagnetism vector I: 

Further, it is convenient to introduce angular variables 
for the vector I. We choose the easiest axis y as a polar 
one, and put 



ly = cost 



Iz = sin 6 cos if , Ix ~ sin 6 sin ip 



Let us set the field in plane perpendicular to the easy 
axis, H II X, then y is always the easiest axis, and there 
is no spin-flop transition. However, at H = He, where 



He - 2{2JDiSY^^ , 



(4) 



the axes x and z switch their roles: the easy plane is xy 
when H < He, and yz when H > He- In this paper, we 
consider the case of strong fields H Ri He; for TMMC the 
critical field is of the order of 100 kOe. In what follows, 
we can safely neglect any contribution coming from D2 
because in strong fields the field-induced anisotropy in 
the basal plane is much larger than the initial crystal- 
field anisotropy. Hereafter, we assume that the field is 
close to He, and regard the quantity 



p = 1 - (H/Hef 



(5) 



as a small parameter playing the role of effective rhom- 
bicity (as far as only static properties are considered, 
H = He corresponds to a "uniaxial" situation). 



III. EFFECTIVE MODEL OF KINK DYNAMICS 

Because of the presence of magnetic field, the equa- 
tions of motion for the model (ph contain not only the 
usual Lorentz-invariant terms, but also terms of the first 
order in time derivatives, and thus their analysis is rather 
complicated. However, one can easily find static soliton 
solutions. There are two types of static kinks, with / ro- 
tating in tbe (xy) plane {ip = ±7r/2) or in the (yz) plane 
{ip — 0, 7r).l3 The behavior of the angle 9 in both kinks is 
given by 



cos 00 = ±tanh(z — z^/zq) 



(6) 



where the kink thickness is zq = Xq for the (yz)-kink and 
zq ~ xoHc/H = xo/{l — pY for the (a;j/)-kink, and Zg 
is an (arbitrary) kink coordinate. For H < He the kink 
with if — ±7r/2 is the fowest energy nonhnear excitation, 
and for H > He this role goes to the kink with 93 = 0, tt; 
ai H = He their energies are equal, as well as the energies 
of in-plane and out-of-plane rnagnon branches. 

As it was shown in Refs. |2^Jl§, approximate soliton so- 
lutions can be constructed using the following trick: One 
can easily see that the characteristic scale of (/? varia- 
tion l^p = xqI sqrtp is much larger than the characteristic 
thickness of the kink x^- Then, "within the kink," i.e., 
within the region where the variation of 9 mainly takes 
place, we can neglect spatial dependence of ip and put 
ip = ^s = const. 

We use the same trick, but regard Zs and Lps as new 
dynamic variables slowly varying in time: we assume that 
the soliton solution has the form 



cos 9g = a tanh 
zq{lp) = xq{1 - 



Zs{t) 



psinV)"'/' 



V = Vs{t) , 



(7) 



here cr = ±1 is the topological charge distinguishing kinks 
and antikinks. Substituting the ansatz ffl) into the La- 
grangian (H) and assuming that Zs and tp^ are "slow" 
variables, i.e. 



^s < c. 



<^s < Wo = c/xo , 



(8) 



one can obtain the following effective Lagrangian with 
only two degrees of freedom: 



L,ff = -i?o(l - psinVs)^/^ + g(i.^ + ^^i.^) 



tt-Eq 
2c 



azg{l ~ p) ' cos Lpg — hSips 



(9) 



Here the dot denotes differentiation with respect to time, 
and Ea = 25*^(2 JI?i)^/^ is a characteristic energy scale 
(rest energy of a soliton a.i H ~ He). The last term in the 
Lagrangian comes from the topological term in (g) and 
is ineffective in classical treatment since it is a complete 
derivative. However, in quantum treatment it changes 
the definition of canonical momentum conjugate to ips 
which is very important, as we will see below. 

The effective potential for the internal variable (ps 
strongly depends on the kink velocity Zg- One can see 
that a classical analysis of the effective Lagrangian (Q) 
allows one to reproduce the results of a more rigorous 
approach of Refs. |l^,|l^,g which we briefly revisit below. 
The canonical momentum P conjugate to Zg, 



E, 



. 



Zs + -Po(l ^ p) 'cr cos yjs , where Pq 



2c 



(10) 



is conserved, and thus it is convenient to express every- 
thing in terms of P. The effective potential 



, TT^Eo r p 



(11) 



8 iPo 



— -acospsil- p) 



1/2 



has four extrema as a function of p>s, at (ps = 0,7r and 
at p>s — ipi^2iP)j where ipi^2iP) are two roots of the 
equation 



coa{ips) = (tP/Pc 



(12) 



Pe^Po{l-p[--- 



1 4 

2 " ^ 



Extrema at ips = 0,tt correspond to (yz)-kinks, and ex- 
trema at ips — ipi.2iP) correspond to {xy)-type kinks. 
When the velocity of an (a;y)-type kink is zero, the plane 
of rotation of the vector I is (xy) , and this plane smoothly 
changes to (yz) when P tends to Pc, so that (xy)-type 
kink smoothly transforms into a (yz)-kink. The fre- 
quency of linear oscillations around those extrema is 
given by 



^ Wo/(p) 



o-Q— - 1 ) for ips =0,^T 



l- 



P2 



(13) 



for (ps = <Pi,2{,P) 



where 



o-Q = (7Sgn[cos(ps) , 



Up) 



T 



1- 



Thus, (a;2/)-type kinks are stable in the entire region 
—Pc < P < Pc where the corresponding solutions ( p^ ) 
have sense. The stability of (j/z)-kinks depends on P and 
on the sign of the quantity ao = sgn(0g cost^s): iyz)- 
kinks with (Jq > are stable at P > Pc, and those 
with tJo < are stable at P < —Pc- For each (j/z)-kink 
there exist two possible "daughter" {xy)-type kinks, cor- 
responding to the two possible solutions (pi^2{P), so that 
the transition from (j/z)-kink to {xy)-]diik when \P\ drops 
below Pc is a spontaneous breaking of symmetry. 

The dispersion relation for soliton excitations, i.e. the 
dependennce of their energy E on the momentum P, is 
schematically shown in Fig. Il| (the curve labeled "clas- 
sical"). For p > 0, i.e. H < He the minimal energy is 
reached at P = 0, while at p < {H > He) the minimum 
is located at P = Pq; in both cases the second derivative 
(PE/dP"^ is discontinuous at P = ±Pc- 

It should be remarked that the effective model (ph was 
derived under the assumption that the variable ips is 
"slow," and therefore expressions (O) have sense only 
if (jjf <^ lOq. This does not matter, however, for the clas- 
sical analysis of stability presented above, because the 
corresponding frequencies are small in the vicinity of the 
transition. But what does matter is (i) strong nonlinear- 
ity: in linear analysis the localized mode frequency is zero 



at P = Pc, which indicates that anharmonic terms have 
to be taken into account; (ii) quantum effects: the magni- 
tude of quantum zero-point fluctuations around the clas- 
sical minima of the effective potential strongly increases 
when the potential well becomes "flat." 



IV. QUANTUM INTERNAL MODES 

Since the effective model (^) contains only two degrees 
of freedom, one can easily quantize it in a canonical way. 
The Lagrangian does not depend on Zg and thus the total 
momentum P is conserved and can be treated as a c- 
number, so that we are left with only one variable ips- 
The Hamiltonian can be easily obtained in the form 



H 



n^u^l 



2Eo { dip. 



d 



S 



Ucsi^s 



(14) 



where S in curly brackets appears because the last term 
in (Q) changes the definition of the angular momentum. 
As we will see below, the presence of this term changes 
the situation drastically depending on whether S is inte- 
ger or half-integer. For the wavefunction ^(</?i) one has 
usual periodic boundary conditions '^{^p) = '^{(p + 27r). 

The effective potential C/off is generally two-well, with 
energetically equivalent well minima; at P = the min- 
ima are located at tps == ±7r/2, as P deviates from zero 
they both move towards or tt depending on the sign 
of P and the 'generalized' topological charge ctq, and at 
P — ±Pc they merge into one well aX ips — Q or tt. One 
may expect that tunneling between equivalent minima 
at \P\ < Pc splits the lowest energy levels, leading to the 
appearance of one more internal soliton mode. However, 
using arguments of the same type as in Refs. |l^j2^, it 
is easy to show that at P = the tunneling is possible 
only for integer S. Indeed, the tunneling amplitude is 
proportional to the path integral 



V2{P) 



■Pl(P) 



Vifis exp < - dtLcff 



(15) 



where (pi^2iP) are the well minima and the effective La- 
grangian is defined by Eq. (^. Since L^s contains the 
topological term hStps, different paths contribute with 
the phase factor proportional to e^^^''^% where Aps is 
the change in ips along the path. At P = one has two 
equivalent paths with Aips = ±7r, so that their contribu- 
tions cancel each other for half-integer S (recall that in 
case of TMMC S" = f). If P 7^ 0, the two paths become 
inequivalent, thus allowing nonzero tunneling amplitude. 
The Schrodinger equation with the Hamiltonian ( |l4| ) 
can be easily solved numerically, and one gets the energy 
eigenvalues as functions of the soliton momentum P; a 
typical level structure is shown in Fig. ^ All states at 
fixed P 



can be labeled with two quantum numbers n and v (y 
being the parity of the "reduced" wavefunction u), and 
for half-integer S odd and even states are degenerate at 
P — Q. One can see that quantum effects change the 
soliton dispersion considerably comparing to the classical 
calculation; particularly, the dispersion tends to have a 
minimum at nonzero P even for fields lower than Hr. 
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FIG. pj. Typical dependence of the soliton energy levels E on 
its momentum P for half-integer spin 5* (levels presented in 
this figure were computed for S = ^); (a) H < He (p = 0.2); 
(b) H > He (p = —0.2). The curve labeled "classical" corre- 
sponds to the result of classical calculatiorjj for the lowest en- 
ergy state; dashed pieces of the curve indicate unstable states. 
For 5* = I the 1 levels are completely inside the magnon 
continuum and thus are also shown with dashed lines. For 
integer S the degeneracy at P = is lifted. 



which means that He is effectively renormalized towards 
its decrease. This renormalization is not small, as seen 
from Fig. |l|, since at p = 0.2 the lowest energy is still 
reached at P 7^ 0. Unfortunately, we are not able to cal- 
culate this renormalization quantitatively: first of all, our 



effective model is valid only at small p, and second, there 
are other quantum effects influencing the effective soli- 
ton energy (the change in the energy of zero-point fluc- 
tuations of magnon modes in presence of a kink) which 
are not considered here; the inteipested reader is referred 
for details to the original papera^jiHj and to the recent 
review .□ The lowest excitations of internal soliton mode 
clearly have purely quantum nature, i.e., they do not 
have any classical counterpart, and are determined by 
the tunnel level splitting. 



V. INTERNAL SOLITON MODES AND 
RESPONSE FUNCTIONS 

In this section we apply the results obtained so far to 
find out how the excitation of internal soliton degrees of 
freedom contributes into the response functions of easy- 
plane antiferromagnet in a strong-field regime (H close 
to He). 

We wish to begin with a remark concerning the relation 
between the gap in the spectrum of linear excitations hwo 
and the kink rest energy Eq in various soliton-bearing 
systems. For the model of structural phase transition 
used by Krumhansl and Schrieffer fiujo oc <i>Oi and Eq cx 
<i>Q, where $o is the magnitude of the order parameter, so 
that near the transition point Eq <C hojQ, which justifies 
neglecting kink-phonon interaction at T <C i?o in Ref. 0. 
For ferromagnets one obtains Eq/Hujo ~ [J/Dy/^, where 
J and D are the exchange and anisotropy constants, re- 
spectively. In case of the Heisenberg magnets J/D ^ 1, 
and Eq ^ hujo, so that the inequality hu>Q '^ T <^ Eq, 
usually assumed in soliton-gas calculations, can be easily 
satisfied. For antiferromagnets the situation is "interme- 
diate": Eq/Hujq ~ 5, and the inequality Eq/Tuvq holds 
only for 5 :^ 1. The relation 5 :^ 1 is often used as 
a condition under which a quantum spin system can be 
treated classically, but in actual practice the strong in- 
equality is not needed, and classical equations of motion 
work well, say, for S = 5/2, as in case of TMMC. Thus, 
in antiferromagnets for realistic S values Eq is only a 
few times greater than Tiloq^ and usually used condition 
TiujQ <^ T <^ Eq is not valid as a strong inequality. (For 
example, in case of strong fields close to He the charac- 
teristic values of kink energy and magnon gap for TMMC 
are i?o — 20 K and hujQ ~ 7 K, respectively). |-||— ■ 

To justify the phenomenological approaclJjl^ to the 
soliton thermodynamics, the only necessary condition is 
smallness of the kink gas density uq comparing to the in- 
verse kink thickness Xq ^ . This means smallness of the ex- 
ponent exp(— i?o/T), i.e., practically, it is sufficient that 
T di ^Eq. Therefore, for antiferromagnets instead of a 
strong inequality Hloq <^ T one must consider the con- 
verse condition Hujq >: T. 

The most general quantity determining the response 
functions is the so-called dynamical structure factor 
(DSF) S°'^{q,uj), which is essentially the space-time 



Fourier component of the two-spin correlation function 
{S"'{n, t)S^{n' , 0)). In case of antiferromagnet, the quan- 
tities of interest are (i) the DSF component at q = 0, 
describing the intensity of response in electron spin res- 
onance (ESR) experiments; (ii) the DSF component at 
q ~ Qb, where Qs is the Bragg wavevector, which is 
usually measured in inelastic neutron scattering (INS) 
experiments. The main contribution to the DSF at 
q = Qb comes from the correlator of the antiferro- 
magnetism vector (Z(z, i)/(z',0), and the component at 
q = is deterniined by the correlator of magnetization 
(m(z, t)rh{z', 0).Lil Soliton contribution to the correlation 
functions can be calculated using the approximation of 
the soliton ideal gasjd which means that at low density of 
kinks the main role is played by one-soliton correlations. 
Then one has to deal with correlators of the form 



dz / dz e 



t iq(z-z') 



X (/(z - z,{t))A{^,,t)f{z' - z,(0)) A(^„ 0)) , (16) 

where A and / are generally certain operators. The 
average should be taken in both quantum-mechanical 
and thermodynamic senses. Quantum-mechanical av- 
eraging is performed with the wavefunctions ^a.p — 
e^^^''ip\{P,(ps) (we use the notation A = {n,i'} for the 
sake of brevity), and the thermodynamic average should 
be performed with the Gibbs distribution fmiction 



w^{P) = (2^?i)-ie-^-(^)/^e-^-(^^^)/^, 



(17) 



where w is normalized to the total density of solitons, 
and the quantity Eg is the change in the free energy of 
the magnon gas due to the presence of a kink, which is at^. 
tributed to the kink in the phenomenological approach:E3 

S,(A,P) =T^ f dkAn^^^{k,\,P)\n[2smh{hul^'^/T)] . 

3 

Here j/ = 1, 2 labels different branches of magnons of the 
continuous spectrum, and lj)^ are the frequencies of cor- 
responding modes. The change in the magnon density 
of states An^^^ (fc. A, P) is connected with the phase shift 

As (fc,A, P) acquired by a magnon of the j-th branch 
with wavevector k after its interaction with the kink hav- 
ing quantum numbers A, P: An^^' = (\/2'K){dl\'s /dk). 
[In fact. Eg contains divergent terms connected with the 
change in the energy of zero-point vibrations between 
the ground state and the one-soliton state. This diver- 
gence, however, is compensated by the counterterms aris- 
ing from the normal ordering of operatorS|rA^ie].ding finite 
quantum correction to the soliton energy.EZTcj] 

Using the above arguments, expression (p^GJ) can be 
rewritten as 

|/(g)P Y. I di^e-(^^(^)+^=(^'^»/^| (A, P|A|A', P~nq)\' 



A, A 



X S{Ex{P) - Ex' [P - hq) - hu) , (18) 



where |A, P) = ■ijj\{P). At low temperatures T <^ Eq one 
may take into account only the contribution of the lowest 
two soliton levels 0^ (see Fig. |l|). 

At r <^ ?iwo, Ss is proportional to exp{—hu;o/T) and 
can be neglected. At T ^ htuo the effect of E^°^* be- 
comes significant; its accurate calculation requires solv- 
ing the soliton-magnon scattering problem to get the 
phase shifts, which is too complicated to be done analyt- 
ically. However, zy-type solitons, i.e. kinks with P > Pc, 
classically have ips — 0,7r, thus being "sine-Gordon" 
in structure, and the phase shifts for the sine-Gordon 
(SG) problem are known. In this case As (k) ~ 
— 2 arctan(fca;o) at small kink velocities f <C c, and for 
T ^ hioo we obtain the well-known SG result 

Ss = -2rin2 - 2T\n{huJo/T) . 

On the other hand, the same SG result should be valid 
for xy-solitons with P — 0. We will assume that the 
same expression can be applied for all P; this assumption 
seems plausible though it cannot be justified in a rigorous 
way. One can expect that possible weak dependence of Ss 
will lead to certain renormalization of the P-dependence 
of the soliton energy levels E\ , but we believe that our 
results will remain qualitatively correct. 



A. Internal modes in ESR 

The static component of magnetization to, according 
to (0), is directed along the Ox axis. Let us assume that 
the pumping rf field is parallel to Oz, then the intensity 
of the ESR signal is proportional to S'^^(0, uj) and is de- 
termined by the correlator of to^ . It is easy to see that 
only the first term in niz, according to (J^, contributes 
to S^^{0,u}) at nonzero frequencies [the last term in (H) 
does not depend on ips and thus yields the delta func- 
tion in UJ, and the contribution from the second term 
just vanishes because at q — the corresponding form- 
factor f{q) in ( p8| ) is zero]. Thus, the significant term 
in niz is proportional to sin Og'ps, and tp^, according to 
(0) is proportional to p^ — hS, where p^p = —ih-^— is 
the operator of the canonical momentum conjugate to 
iPs- Again, hS yields (5(w), so that the only contribution 
at nonzero uj comes from p^p. Typical plots of the ESR 
response depending on frequency for 5* = | are shown 
in Fig. g; one can see that the transitions between 0+ 
and 0~ soliton levels lead to a pronounced peak with a 
maximum located well below the antiferromagnetic res- 
onance frequency luq. Since the levels 0+ and 0^ have 
different minima in their P dependencies, the position 
of peaks in the ESR response at positive and negative uj 
is also different. When the temperature increases, peaks 
shift towards ujq, but still keep their shape with a clearly 
pronounced maximum. Thus we expect that this effect 
can be really observable experimentally in TMMC, where 
at H = He the three-dimensional ordering temperature 
Tjv ~ 3K, i.e. about O.lbEo. 



B. Internal modes in INS 

In neutron experiments, it is usually possible to mea- 
sure either longitudinal (with respect to the magnetic 
field) DSF component S" = S^^{q,uj), or the transverse 
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FIG. 0. Plots of S^^ {0,<-^) (this quantity determines the inten- 
sity of the ESR signal) for S = | and different values of the 
parameter p=l~ H^/H^: (a) T = 0.05-Bo; (b) T = O.lS^o- 
Arrows show the positions of the boundary of the continuous 
magnon spectrum at _ff = He- 
one S^ = Syy{q,uj) + S''^{q,uj), with q close to the Bragg 
wavevector Qb- In TMMC, both those components ex- 
hibit presence of the soliton central peak.lia In the trans- 
verse DSF the dominating contribution comes from S^^ 
(which corresponds to the contribution of flips of the di- 
rection of the antiferromagnetism vector I along the easy 
axis Oy) , and S'^'^ is insensitive to the effects of internal 
modes because ly does not depend on ips] that is why 
we will be interested only in the longitudinal component 
S^^ determined by the correlator of Ix ■ Fig. H shows sev- 
eral typical plots of S'^'^ {Q b , ijj) for different values of the 



field and temperature (actually, we have also calculated 
S^^{q,u}) for a few values of q close to Qs, but the re- 
sulting curves exhibit very weak q dependence, so that we 
restrict ourselves to presenting only plots for q = Qb)- 
One can see that the most significant manifestation of 
internal soliton modes is a nonzero width of the central 
peak (CP) 
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dominating central peak (see Fig. ^d) , and in TMMC at 
H Ki He ioi T >: T/v — O.lSE'o this peak should be almost 
completely suppressed. 



VI. SUMMARY 

We consider the problem of soliton dynamics in Id 
"almost-classical" spin-S* easy-plane antiferromagnet in 
strong external magnetic field, taking into account in- 
ternal degrees of freedom of solitons. We show that 
quantum effects of tunneling between two energetically 
equivalent kink states lead to level splitting for all values 
of S; the usual selection rule which prohibits the tun- 
nelling for half-integer S is lifted due to the strong cou- 
pling between internal and translational modes. It is pre- 
dicted that this internal mode can be detected in TMMC 
[(CD3)4NMnCl3] by means of the electron spin resonance 
(ESR) or inelastic neutron scattering (INS) technique; we 
also show that ESR is a preferable method because the 
corresponding resonance peak is much more pronounced 
comparing to that in the INS response. 
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FIG. g. Plots of the longitudinal DSF S" 
and different values of the parameter p 
T^OmEo; (b) T = 0.1Eo. 



'{Qb,uj) for 5: 
= 1 - H'^/H'^: 



at g = Qb (recall that in the standard ideal-gas 
calculationQ without taking into account internal degrees 
of freedom the CP width for the longitudinal component 
is proportional to \q — Qb\)- However, it should be men- 
tioned that the effects of soliton-magnon and soliton- 
soliton scattering, which we do not consider here, also 
lead to nonzero CP width at g = Q_bJ23 and therefore 
it may be problematic to separate those two contribu- 
tions experimentally. The additional peak at nonzero 
frequency is also present in 5^^, but with increasing 
temperature it quickly gets smeared by the tails of the 
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FIG. 1. Typical dependence of the soliton energy levels 
E on its momentum P for half-integer spin S (levels pre- 
sented in this figure were computed for S — ^): {a) H < He 
(p = 0.2); (b) H > Hc{p= -0.2). The curve labeled, "classi- 
cal" corresponds to the result of classical calculationtl for the 
lowest energy state; dashed pieces of the curve indicate un- 
stable states. For 5=1 the 1* levels are completely inside 
the magnon continuum and thus are also shown with dashed 
lines. For integer 5* the degeneracy at P = is lifted. 



FIG. 2. Plots of S^^ (0, ix)) (this quantity determines the in- 
tensity of the ESR signal) for S = | and different values of the 
parameter p=l- H^/Hl: (a) T = 0.05.Bo; (b) T = 0.15^0. 
Arrows show the positions of the boundary of the continuous 
magnon spectrum ai H = He. 



FIG. 3. Plots of the longitudinal DSF S'^="(Qs,c^) for 
S — ^ and different values of the parameter p = 1 ~ H /H^ : 



(a) T = 0.05£o; (b) T = O.l^o. 



